Influence of solvent granularity on the eff'ective interaction between charged colloidal 

suspensions 
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We study the effect of solvent granularity on the effective force between two charged colloidal 
particles by computer simulations of the primitive model of strongly asymmetric electrolytes with 
an explicitly added hard sphere solvent. Apart from molecular oscillating forces for nearly touching 
colloids which arise from solvent and counterion layering, the counterions are attracted towards 
the colloidal surfaces by solvent depletion providing a simple statistical description of hydration. 
This, in turn, has an important influence on the effective forces for larger distances which are 
considerably reduced as compared to the prediction based on the primitive model. When these forces 
are repulsive, the long-distance behaviour can be described by an effective Yukawa pair potential 
with a solvent-renormalized charge. As a function of colloidal volume fraction and added salt 
concentration, this solvent-renormalized charge behaves qualitatively similar to that obtained via 
the Poisson-Boltzmann cell model but there are quantitative differences. For divalent counterions 
and nano-sized colloids, on the other hand, the hydration may lead to overscreened colloids with 
mutual attraction while the primitive model yields repulsive forces. All these new effects can be 
accounted for through a solvent-averaged primitive model (SPM) which is obtained from the full 
model by integrating out the solvent degrees of freedom. The SPM was used to access larger colloidal 
particles without simulating the solvent explicitly. 



PACS: 82.70.Dd, 61.20.Ja 
I. INTRODUCTION 

Most of soft matter systems, as colloids, polymers or 
biological macromolecules, are dispersed in a molecu- 
lar solvent Q. Therefore, a full statistical description 
of supramolecular solutions should include the solvent 
explicitly. Such a treatment is highly non-trivial, how- 
ever, since the length scale separation between the meso- 
scopic particles and the molecular solvent directly implies 
that the number of solvent particles which have to be in- 
cluded is enormous. On the other hand, one is interested 
mainly in properties of the big particles such that a sol- 
vent pre-average makes sense. The crudest form of such a 
course-grained level is to treat solvent properties just by 
a dielectric background or by some effective parameters 
which enter in the effective colloidal interactions. This 
procedure is questionable for polyelectrolytes where the 
solvent couples directly to the counterions which may af- 
fect the effective interaction between the polyelectrolytes 
via the long-ranged Coulomb coupling of counterions to 
the polyelectrolytes. 

In this paper, we consider the case of two spherical 
charged colloidal particles (polyions) which are immersed 
in a bath of the molecular solvent and their oppositely 
charged counterions plus additional salt ions Our 
main focus is the total effective force acting onto the col- 
loidal pair which is the key quantity to understand col- 
loidal stability and which governs colloidal correlations 
and phase transitions. In almost any theoretical treat- 
ment, the discrete structure of the solvent particles was 
neglected and only the charged constituents were treated 
explicitly within the so-called "primitive" model (PM) of 



strongly asymmetric electrolytes. Even this model is non- 
trivial in the colloidal context due to the large asymmetry 
between poly- and counterions and bears a rich physics 
resulting from the strong coupling between the different 
species. In recent computer simulations j^-Ql, counteri- 
onic correlations have been shown to be responsible for ef- 
fective attractions between the like-charge polyions. The 
PM, re-formulated in terms of modern density functional 
theory of the inhomogeneous counterion plasma can 
also also used as a starting point to derive simpler theo- 
ries such as the mean-field nonlinear Poisson-Boltzmann 
approach or the linearized Debye-Hiickel-type screening 
theory. The latter results in an effective Yukawa pair 
potential between the colloids as given by the electro- 
static part of the celebrated Derjaguin-Landau-Verwey- 
Overbeek (DLVO) theory [^. This potential can also 
be used with renormalized parameters to include parts 
of the nonlinear screening effects arising from Poisson- 
Boltzmann theory |]l0| . 

In the present paper, we investigate the influence of 
solvent granularity on the effective interactions between 
charged colloids. We model the solvent as a hard sphere 
fluid at intermediate packing fractions and use computer 
simulations and the theoretical concept of effective in- 
teractions to derive effects due to the discrete solvent. 
The PM is tested against this more general model. Al- 
though the hard sphere model neglects some important 
solvent properties as its polarizibility and its per- 
manent multipole moments |12|, it provides a minimal 
framework to get insight into counterion hydration and 
screening effects. The hard sphere solvent model (which 
is sometimes called solvent-primitive model) has been 
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used also in many other investigations of ordinary elec- 
trolytes and for electrolytes confined between two par- 
allel charged plates. Most of the approaches invoke ad- 
ditional approximations as different versions of liquid- 
integral equations Jl3| , p^ , p^ , Poisson-Boltzmann theory 
suitably modified to include the short-ranged solvent de- 
pletion effects or more sophisticated density func- 
tional approaches of multicomponent systems |lq ]. For 
charged plates and for small neutral particles |l^] 
some computer simulations have already been performed 
including a hard sphere solvent explicitly but there are 
no results for charged colloidal spheres. 

Most of the results in this paper are based on a new 
"solvent bath" simulation scheme which allows to sim- 
ulate many neutral spheres together with the charged 
species. For divalent counterions, we obtain attractive 
forces due to overscreening of polyions by counterions 
which are attracted towards the colloidal surfaces via 
hydration (or solvent depletion) forces. For monovalent 
counterions and large distances we show that the con- 
cept of charge renormalization can be used to extract a 
Yukawa picture of the effective interaction with a solvent- 
renormalized polyion charge. We check the trends of this 
renormalized charge with respect to the colloidal den- 
sity and the concentration of added salt and find qualita- 
tive agreement but quantitative differences as compared 
to the Poisson-Boltzmann theory. All our results can 
be reproduced within a solvent-averaged primitive model 
(SPM) which was extensively used in earlier theoretical 
studies of electrolytes between plates JTs],^ . This idea 
originates from McMillan and Mayer pO(| dating back to 
1945. 

Our paper is organized as follows: In chapter II, we 
describe our model and define approximations on differ- 
ent levels. The computer simulation method is described 
in chapter III. We then turn to results for the neutral 
case in chapter IV and for the salt-free case in chapter 
V. Parts of the latter have been published elsewhere . 
The effect of added salt is described in chapter VI and 
other mechanisms of polyion-polyion attraction are crit- 
ically discussed in chapter VII. We finally conclude in 
chapter VIII. 



II. MODELING ON DIFFERENT LEVELS 

In this section we summarize the modeling on different 
levels. In the following we shall use the most detailed 
description of the hard sphere solvent model and test the 
validity of the different inferior levels with respect to our 
data. 



A. The hard sphere solvent model (HSSM) 

The hard sphere solvent model (HSSM) involves 
spherical polyions with diameter <Tp and homogeneously 



smeared charge qp together with their counterions of di- 
ameter (Tc and charge qc in a bath of a neutral solvent 
{qs — 0) with diameter as- In the absence of salt, the 
pair potentials between the particles as a function of their 
mutual distances r are a combination of excluded volume 
and Coulomb terms 



oo 



for r < (cTi 
else 



(1) 



where e is the a smeared background dielectric constant 
of the solvent and {ij) = {pp), ipc), (ps), (cc), (cs), (ss). 
Further parameters are the thermal energy ksT and the 
partial number densities pi {i — p, c, s) which can be 
expressed as partial volume fractions (j>i = irpiaf/G {i = 
p,c,s). Charge neutrality requires = pd^cl- Addi- 

tional salt ions can readily be included into the descrip- 
tion as further charged hard spheres. 



B. The solvent-averaged primitive model (SPM) 

For a fixed configuration of charged particles the sol- 
vent can be traced out exactly arriving at depletion forces 
F^'^'' acting onto the ith charged particle. They can be 
related to a surface integral over the solvent equilibrium 
density field Psif) which depends parametrically on the 
positions of the fixed charged particles: 



^ksT I df ps{r^ 

Si 



(2) 



where / is a surface vector pointing towards the center 
of the ith charged particle. If one adds these forces to 
the PM, the resulting model is strictly equivalent to the 
HSSM. The integrand Ps{r) is affected by the space ex- 
cluded for the solvent due to the presence of the finite 
core of the charged particles resulting in inhomogeneous 
density distributions around the excluded volume. The 
range of this inhomogeneity is characterized by the hard 
sphere bulk correlation length ^ which depends on the 
solvent packing fraction 0s . A further approximation de- 
composes the forces Pj^'^'^ into pairwise parts, i.e. into a 
superposition of pair contributions coming from neigh- 
boring charged particles. This approximation is justi- 
fied if the average distance between triplets, quadruplets, 
etc. of charged particles is much larger than the bulk 
correlation length ^. In the salt-free case, this is gen- 
erally granted except for nearly touching polyions with 
squeezed counterions. If salt is added, the justification 
is less clear as ion pairing by counter- and colons near 
polyions may be an important configuration. 

The resulting pairwise interactions define the solvent- 
averaged model (SPM) where the depletion pair poten- 
tials V^j'^\r) {{ij) — (pp), {pc), (cc)) have to be added to 
the interactions of the primitive model of the next para- 
graph. These pairwise depletion forces have been the 
subject of intense recent research |22-27|. In particular, 
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we will determine them by computer simulation, and use 
these results as an input for the SPM. 



C. The primitive model (PM) 

The primitive model has the same interactions as the 
HSSM except for the absence of the solvent. Hence the 
basic interactions are again 

Too forr<(a.+a,)/2 

I QiQj/'^^ else ^ ^ 

but now for (ij) ~ (pp), (pc), (cc) only. 



D. DLVO-theory 

In DLVO theory, only the polyions are treated explic- 
itly. The electrostatic part of their interaction is an ef- 
fective Yukawa pair potential which has the form: 



Vir) 



gp exp(-K(r - o-p)) 
(1 + Kcrp/2)2er 



with 



ft = yj4,TTpcqllekBT 



(4) 



(5) 



E. The PB-renormalized Yukawa model (PBYM) 

This approach was suggested by Alexander et al | |Tc| ] 
and is based on Poisson-Boltzmann theory in a spherical 
cell around a single polyion. The cell radius R is fixed 
by the polyion concentration: 



R - (47rpp/3)-i/3 



(6) 



Within Poisson-Boltzmann theory, one calculates the 
countcrion density pc at the cell boundary. Lineariz- 
ing the nonlinear Poisson-Boltzmann theory at the cell 
boundary, one obtains again an effective Yukawa poten- 
tial between the colloids arriving at the PB-renormalized 
Yukawa model (PBYM). The Yukawa potential has the 
same form as in Eqn.(H) but contains a rcnormalizcd in- 
verse screening length 



and a renormalized charge 



(7) 



(8) 



which is considerably smaller than the bare charge qp. 
Many experimental data for the colloidal structural cor- 
relations 123] , the long-time self-diffusion or the freez- 
ing line pOf] have been analyzed using this concept of 
charge renormalization and in general good agreement 
was found for monovalent counterions provided the col- 
loids are far away from charged plates |3l[|. 



F. The solvent-renormalized Yukawa model (SYM) 

This approach is a generalization of the Poisson- 
Boltzmann cell modcll ijl^l to the presence of a gran- 
ular solvent. Again, one considers a single polyion 
in a spherical cell, but uses the full HSSM to obtain 
the countcrion density pc at the cell boundary. As in 
the PBYM, the associated Yukawa pair potential has 
a solvent-renormalized inverse screening length k* and 
a solvent-renormalized charge g*, which, however, differ 
from that of the PBYM approach. 



III. SIMULATION METHOD 

We consider two large spherical polyions in a cubic 
simulation box of length L with periodic boundary con- 
ditions, hence pp = The polyions are fixed along 
the body diagonal of the cubic box. While the simu- 
lation methods for the PM are straightforward and are 
described elsewhere [^H), a significant volume fraction 
of solvent particles together with large colloidal particles 
implies a huge number of solvent spheres in the simula- 
tion box. 

As the solvent interactions are short ranged and the 
solvent-averaged interactions are only of range ^, one 
can reduce the number of solvent particles in the sim- 
ulation box considerably using a "solvent-bath" method. 
This procedure is sketched in Figure |] and works as fol- 
lows: we define a spherocylindrical cell around the col- 
loidal pair such that the minimal distance h from the col- 
loidal surface to the spherocylindrical boundary is much 
larger than the hard sphere bulk correlation length ^. 
The hard sphere solvent is only contained in this sphero- 
cylinder. As the cell volume is considerably smaller than 
the volume of the whole simulation box, the num- 

ber of solvent particles for fixed given volume fraction 
(f>s is drastically reduced. No restriction is done for the 
counter- and salt-ions, which can move within the whole 
simulation cell. 

Care has to be taken at the artificial cell boundary. 
We use the Molecular Dynamics (MD) method calcu- 
lating the particle trajectories and performing statistical 
averages over some physical quantities. The hard sphere 
solvent is treated by the well-known hard sphere collision 
rules. Once a solvent particle is leaving the spherocylin- 
dric cell it is randomly inserted at another place of the 
cell boundary with the same velocity. The position of its 
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random insertion has a randomly chosen distance from 
the spherocyhndrical boundary up to which avoids 
unphysical solvent layering there. 




FIG. 1. View of the set-up as a downward projection of 
a simulation snapshot: Two polyions (dark open circles) in 
a bath of solvent particles (small hollow spheres) contained 
in a spherocylindric cell of width h. The counterions shown 
as small dark spheres can move in the whole simulation box 
of size L. The distance between the polyions is r, hence the 
projected distance shown is r/%/3. 

Since the width of the cell h is much larger than the 
hard sphere bulk correlation length ^, the presence of 
the boundary has no influence on the inhomogeneous 
density distribution of the solvent and the counterions 
near the colloidal surfaces. The counterion motion is im- 
plemented as follows: Outside the spherocylindrical cell 
the counterions interact via a solvent-averaged effective 
depletion potential. This is justified as the typical dis- 
tance between the counterions is much larger than the 
hard sphere bulk correlation length ^. Therefore, the 
correction to Vccir) due to solvent layering is negligibly 
small anyway for counterions outside the spherocylindri- 
cal cell. This is not the case when salt is added as the 
attraction between co- and counterions may lead to short 
distances where solvent depletion effects may dominate 
the interactions. For a counterion approaching the cell 
boundary, there is an artificial asymmetry between the 
solvent bath inside the cell and the "vacuum" outside the 
cell which hinders a counterion to penetrate into the sol- 
vent bath. This unphysical effect is repaired in the sim- 
ulation scheme by switching off the counterion-solvent 
interaction for a counterion which is penetrating from 
outside. Once the counterion is fully surrounded by sol- 
vent molecules the interaction is turned on again. This 
leads to a symmetric counterion crossing rate across the 
spherocylindrical cell boundary. 



We have carefully tested the algorithm against simu- 
lations of small systems where the whole space was filled 
with solvent particles. We also tested against a situa- 
tion of a solvent slab between charged plates. Perfect 
agreement was found compared to simulations where the 
whole space was filled with solvent particles. 



IV. RESULTS FOR THE NEUTRAL CASE 

Let us first discuss the much simpler case of neutral 
polyions {q^ — 0) under absence of counterions. The re- 
sulting system is just a pair of big hard spheres in a sea 
of small solvent spheres. This simple model for binary 
mixtures of hard sphere colloids has gained considerable 
attention during the past ten years. The effective inter- 
action between the large spheres as induced by depletion 
of the small spheres in the zone intermediate between 
the nearly touching big spheres has the following char- 
acteristic features which were obtained by density func- 
tional theory |p3|~p5|| , computer simulation ||3^ , |26| , p7|| and 
experiments pq-pW: it is attractive for nearly touching 
spheres, then it oscillates with the bulk correlation length 
of the hard sphere solvent ^. The depletion interactions 
decays to zero exponentially with the surface-to-surface 
separation of the big spheres. The characteristic decay 
length is again the bulk correlation length ^. 

Our motivation to investigate the neutral case is 
twofold: First, it is a simple case which allows to test 
our simulation set-up. There are some computer simula- 
tions of the depletion interaction in the literature but the 
range of parameters examined is far from being complete. 
Second, for the solvent-averaged primitive model (SPM) 
it is exactly the solvent depletion term which one has 
to add to the primitive interactions as given by eqn.(|^). 
Therefore, studies of the SPM require a full knowledge 
of the neutral case. 

Computer simulation results for the total effective de- 
pletion force Fpp-* (r) acting onto the big spheres are pre- 
sented in Figure for two size ratios (Jp/us of 2 and 14. 
The solvent packing fraction was chosen to be (j)s = 0.3. 
The force is projected onto the particle separation vec- 
tor such that a positive sign means repulsion. We note 
that a direct evaluation of eqn. (||) is difficult as the 
solvent density field piles up strongly at the surfaces of 
the big particles. A much more effective way is to mea- 
sure the momentum transfer on the fixed big spheres due 
to colliding small spheres during the Molecular Dynam- 
ics simulation. The results are compared with a recently 
proposed fitting formula of Roth and Evans based on 
density functional theory . One sees that the contact 
value of the force and the oscillations are well-described 
by the theory but there are deviations around the first 
maximum. This is less apparent if the effective potential 
is compared as shown in the inset of Figure ^. Obviously, 
the reason for that is that the force is a derivative which 
is more sensitive to approximations. In the computer 
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FIG. 2. Reduced distance-resolved depletion force F{r)as/kBT versus reduced distance r/us between two identical neutral 
spheres embedded into solvent bath of packing fraction (j)s = 0.3: a) size asymmetry of cTp : cts = 2 : 1; b) size asymmetry of 
(jp : (Js = 14 : 1. Solid line- our simulation results, dashed line- fitting formula of Ref. [38]. The inset shows the corresponding 
reduced depletion potentials Vp'p {r) /ksT . 



simulation, the effective potential Vpp\r) can be accessed 
by integrating the distance-resolved computer simulation 
results of the force 



dr'F^tHr') 



(9) 



Also the whole set of depletion pair potentials V^^^ (r) 
{{ij) = (pp) , (pc) , (cc)) which are the input of a typ- 
ical SPM simulation are presented in Figure |. The 
counterion-counterion interaction is dominated by the 
Coulomb repulsion also shown as a dashed line in Figure 
3c. The bare Coulomb repulsion between the polyions 
is much larger than the polyion-polyion depletion poten- 
tial and is not shown in Figure 3a. Finally, the polyion- 
counterion depletion interaction exhibits a deep attrac- 
tion near contact of the order of ksT which is of sim- 
ilar order than the Coulomb attraction also shown as 
a dashed line on Figure 3b. This will have important 
consequences of counterion adsorption on the colloidal 
surface. This effect is induced by the granularity of the 
solvent and is absent in the PM. 



V. RESULTS FOR THE SALT-FREE CASE 

A. Nano-sized colloids 

Although the amount of solvent which has to be sim- 
ulated explicitly has been reduced drastically by the 
solvent-bath scheme, only colloidal sizes which are in the 
nano-domain can be addressed on present-day comput- 
ers. We have performed extensive computer simulation in 



this domain to check carefully the different approaches. 
We find that the SPM describes the full simulation data 
of the HSSM very well. Larger colloidal sizes are thus 
only accessible within the SPM and discussed in chapter 
V.B. 

In our simulations, we fixed T — 298°K and e = 81 
(water at room temperature) with as — 3A, 0s = 0.3 
such that f is about Sag- We varied the polyion charge 
and size and the counterion diameter CTc. The width of 
the spherocylindrical cell h is IOcTc such that typically 
Ng — 25.000 — 30.000 solvent hard spheres are simulated. 

We have basically calculate two quantities: first, as 
a reference, we have calculated the spherically averaged 
counterion density profile Pc{r) around a single polyion 
where r is the distance from the polyion center. The sim- 
ulation was done in a cubic box of reduced length 
with periodic boundary conditions in order to reproduce 
the colloidal packing fraction (j)p. Second, our target 
quantity is the total force F{r) acting onto a polyion for 
a given colloid-colloid separation r. This effective force 
F(r) is the sum of four different contributions: 

i) the direct Coulomb repulsion as embodied in Vpp{r) 
(note that all the periodic images contribute to the total 
force) , 

ii) the counterion screening resulting from the averaged 
Coulomb force of counterions acting onto the polyions, 

iii) the counterion depletion term arising from the hard 
sphere part of Vpc{r), 

iv) the solvent depletion force. 

Explicit results for F{r) are presented in Figure ^ 
where the solvent and the counterion diameter were cho- 
sen to be equal and the counterion were monovalent. 
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FIG. 3. Reduced depletion potentials V^f\r)/kBT 
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= (pc), (cc)) versus reduced distance r/os: 

a) polyion-polyion depletion with a^jos ~ 10; b) 
polyion-counterion depletion with ap : ac o^s ~ W : 1 : 1 
a) counterion-counterion depletion with ad a a = 1 The sol- 
vent packing fraction is (^s = 0.3. The inset shows the sit- 
uation. The dark spheres correspond to the pair of charged 
particles, the solvent are the hollow spheres. The a;-axis starts 
for touching particles. For comparison we have also included 
the Coulomb interaction of the PM as dashed lines in b) 
and c) for = — 32e,gc = le and e = 81. Note that, the 
polyion-counterion Coulomb potential in plot 3b is reduced 
by a factor of 1/10. 




FIG. 4. Reduced distance-resolved force F(r)(Tp/kBT ver- 
sus reduced distance r/op. The inset shows the same for 
nearly touching polyions of molecular distances. The sim- 
ulation parameters are: qc ~ le, qp = — 32e, e = 81, 
(Tp : (Tc : (Ts = 10 : 1 : 1, = 5.8 X 10"'^. Solid line with error 
bars: -HSSM; long-dashed line: SPM; short-dashed line: PM; 
open circles: DLVO theory; open squares: PBYM theory; 
dotted line : SYM theory; dot-dashed line in inset: solvent 
depletion force (for comparison). 

The force exhibits oscillations for molecular distances 
due to solvent and counterion layering and is repulsive 
for larger distances. The SPM yields surprising agree- 
ment with the HSSM describing even the molecular os- 
cillations for nearly touching polyions, sec the inset of 
Figure | , while the PM overestimates the force consider- 
ably. This can be attributed to the fact that the SPM 
incorporates the additional counterion accumulation at 
the colloidal surface due to the hydration or solvent de- 
pletion. This can clearly be seen in the counterionic den- 
sity profile around a polyion as shown in Figure ^ which 
piles up near the colloidal surface. While this accumula- 
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tion is quantitatively described by the SPM it is absent 
in the ordinary PM. The PBYM and DLVO theory lead 
to forces which strongly overestimate the HSSM data. 
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FIG. 5. Reduced counterion density profile pcO^ around a 
single polyion versus reduced distance r/up from the polyion 
center. The parameters and the line types are as in Figure ^ 

We have further tested the frequently invoked ^''super- 
position principle^'' which approximates the total force as 
a sum of the PM and the depletion term. Its comparison 
to the full HSSM data is given in Figure 6a. The first 
maximum of the total force is semi-quantitatively repro- 
duced but the superposition principle predicts a second 
maximum which is too sharp as compared to the HSSM 
data. This becomes even worse for a doubled counte- 
rion diameter of 6A where the superposition predicts a 
secondary maximum which is completely absent in the 
HSSM data, see Figure 6b. The physical reason for that 
is that the counterion layering coupled to the solvent de- 
grees of freedom becomes relevant for these distances. 

The forces for a doubled counterion diameter Uc are 
presented in Figure 0. For small distances (except for 
touching), the PM yields larger forces as compared to 
Figure 4, as the counterion repulsion is stronger which 
reduces screening. In the HSSM and SPM, on the other 
hand, also the polyion-counterion depletion attraction is 
getting stronger, such that the total polyion screening is 
practically unaffected. Of course, the PBYM and DLVO 
theory yield results which are insensitive to the counte- 
rion diameter. 

Furthermore we have investigated the case of stronger 
Coulomb coupling by considering divalent counterions. 
Explicit data are shown in Figure There is overscreen- 
ing of polyions resulting in a mutual attraction between 
like-charged polyions. We emphasize that it is the elec- 
trostatic term of the counterions that produces the at- 
traction but not the counterion or solvent depletion term. 
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FIG. 6. Test of the superposition principle: reduced dis- 
tance-resolved force F{r)ap/kBT versus reduced distance 
r/(7p for two different counterion sizes: a) CTc — 3A, b) 
(7c — 6A. The force predicted by the superposition principle 
is the short-dashed line. The other parameters and notations 
are the same as in Figure ^. 

Nearly every counterion is in the presence of the colloidal 
surfaces, as demonstrated by the countcrionic density 
profile shown in Figure ^ where the piling-up of coun- 
terions near the colloidal surface is much stronger. The 
attractive force has a range of several polyion diameters. 
Again the SPM perfectly reproduces the forces. The PM 
(and also the PBYM and DLVO theory), on the other 
hand, yield repulsion. This demonstrates that a discrete 
solvent has a profound infiuence on the effective interac- 
tions. 
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FIG. 7. Same as in Figure ^ but now for a double counte- 
rion diameter such ap : a"c : cts = 10 ; 2 : 1. 




1.3 2 3 4 5 

FIG. 8. Same as in Figure ^ but now for divalent counteri- 
ons and ap : ac <Js = 14 : 2 : 1. The further parameters are 
\qp/qc\ = 32 and (t>p = 5.8 x 10"^ 

We finally discuss the validity of the solvent- 
renormalized Yukawa model (SYM). Computer simula- 
tions have been performed for a single polyion in a spher- 
ical cell and the counterion boundary density was calcu- 
lated. The boundary of the cell was not hard but coun- 



terions leaving the cell were inserted at the opposite side 
of the cell. Again a smaller spherical solvent bath around 
the polyions with a width h was used, see Figure |o| for 
the set-up and a projected simulation snapshot. As Fig- 
ures 4 and 7 show, the SYM is indeed a reasonable de- 
scription of the forces for large distances. We also remark 
that the SPM and the HSSM yield the same counterion 
density at the boundary of the spherical cell needed as an 
input for the SYM which justifies the usage of the SPM to 
get the solvent-renormalized Yukawa parameters of the 
SYM. 

1 n , . , . , . , 1 




L ^ ' 

0.6 0.65 0.7 0.75 




FIG. 9. Same as Figure 5 but now with the parameters of 
Figure || 

The validity of the SYM only holds for the case of mono- 
valent counterions where the remaining "free" countcri- 
ons are responsible for the screened, repulsive force. For 
divalent counterions, no free counterions are left, and a 
linearized screening theory breaks down such that an at- 
traction cannot be encaptured by the SYM. 

B. Mesoscopically-sized colloids 

Our results for mesoscopically-sized colloids are based 
on SPM simulations as justified in the previous chapter. 
Distance-resolved colloidal forces F{r) for monovalent 
counterions, a size asymmetry of Cp : CTc : CTs = 370 : 1 : 1 
or 370 : 2 : 1, and a charge ratio of Qp/qc = 280 are 
presented in Figure These forces are repulsive but 
much smaller than that from PM simulations. Again, 
this is due to counterion accumulation near the colloidal 
surface as induced by the additional solvent depletion at- 
traction. As the corresponding potential energy gain is 
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only few fc^T, this depletion attraction is different from 
chemisorption of counterions. The solvent-renormalized 
Yukawa model (SYM) leads to forces which are very simi- 
lar to the SPM over the whole range of distances explored 
while the PM overestimates the forces. 




FIG. 10. View of the set-up and projected simulation snap- 
shot for a single polyion in a spherical Wigner-Seitz-cell. The 
polyion is shown as dark open circle in the cell in a bath of 
solvent particles (small hollow spheres) contained in a spher- 
ical cell of width h. The counterions are shown as small dark 
spheres. R is cell radius. 

The traditional meaning of the "bare" charge qp in the 
PM is not the full polyion charge but a smaller charge 
which results from a polyion charge reduction by strongly 
adsorbed (or condensed) counterions. This picture can 
also be tested against our results. We first have calcu- 
lated the average number of counterions in a molecu- 
lar shell around the colloids of width ^. If the polyion 
charge is reduced by this amount and the PM is used to 
predict the effective interaction, the resulting force still 
overestimates the HSSM data, see the open diamonds in 
Figure 11a. In order to fit these data satisfactorily, one 
has to assume an unphysically large width of 5^ to get a 
charge reduction that reproduces the SPM data. Hence 
the PM cannot be justified even with a polyion charge 
reduction. The reason for that are the weak hydration 
forces which are quite different from chemisorption pro- 
viding a strong counterion binding with an energy gain 
of hundreds or thousands of ksT. Furthermore, an ar- 
bitrary splitting into a fraction of condensed counteri- 
ons and "free" counterions described by DLVO, Poisson 
Boltzmann or any other local density functional theory is 
not possible: near the colloidal surface the electric dou- 
ble layer is highly correlated such that fixing a fraction of 
counterions gives a completely different picture. Only if 
the fraction of free counterions is determined within an 



approach that includes all these correlations (as in the 
SYM), a linearized screening theory far away from the 
colloidal surfaces is justified. 



400 rr 




12 3 4 




FIG. 11. Reduced distance-resolved force -F(r)ap/fc_BT ver- 
sus reduced distance ^/<^p 
for larger polyions, ap : ac : cfs = 370 : 1 : 1, \qp/qc\ = 280, 
cj>p — 2.3 X 10""^, and monovalent counterions. a) Long-dashed 
line: SPM; short-dashed line: PM; open circles: DLVO the- 
ory; open squares: PBYM; open diamonds: PM with charge 
reduction, dotted hne : SYM. b) Long-dashed line: SPM for 
a doubled counterion diameter ac = 6A; solid line: SPM for 
ac — 3A; dashed line: PM. The inset shows the corresponding 
reduced counterion density profile in the vicinity of a single 
polyion versus reduced distance d/as, d being the distance 
from polyion surface. 
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This consideration casts some doubts to some recent 
theories where such a partitioning of counterions is an 
input, see e.g. Refs. [|9|,|o). Further SPM results for a 
doubled counterion diameter are presented in Figure lib. 
As the counterion depletion force is getting stronger for 
a large counterion, the force is getting smaller, compare 
the full and dashed line in Figure lib. The PM (short- 
dashed line in Figure lib), on the other hand, is practi- 
cally insensitive to a change of the counterion diameter 
except very close to the colloidal surfaces. This picture 
gains further support from the counterionic density pro- 
files around a single polyion shown in the inset of Figure 
lib for distances very close to the colloidal surface. A 
layer of condensed but still mobile counterions close to 
the surfaces is present in the SPM which is absent in 
the PM. The larger the counterion diameter the more 
counterions are in this layers as the depletion is getting 
stronger. 




150 



100 I ' ' ' ' ' ' ' 1 

-4 -3 -2 -1 

FIG. 12. Renormalized charge versus decadic logarithm 
logjQ (j)p of the polyion fraction as obtained within in spherical 
cell containing a single polyion. The parameters are the same 
as in Figure |ll[ Solid line: SPM; long dashed line: PBYM; 
dashed line: PBYM for a fixed bare charge of qp — 269e. 

We finally discuss the solvent-renormalized charge q* 
as a function of the colloid volume fraction (j)c for fixed 
bare charge qp and compare with the prediction of the 
traditional charge renormalization approach within in 
Poisson Boltzmann cell theory (PBYM) |l0) . Simulation 
data for q* based on the SPM in a spherical cell are shown 
on the full line in Figure The renormalized charge is 
smaller than the bare charge and behaves non-monotonic 
with the particle density. The non-monotonicity is sta- 
ble with respect to added salt and is related to a non- 
monotonic counterion density at the cell boundary as a 
function of density. It can be understood as follows: For 
extremely high packing fractions the spherical cell acces- 



sible for the counterions is a very thin shell across which 
the polyion-counterion attraction varies slowly. Due to 
the rapidly decreasing volume accessible for the counteri- 
ons, the boundary counterion density becomes larger for 
increasing (/)p, see the volume fraction correction in Refs. 

On the other hand, for very small (j)p, entropy 
of counterions will force them to cover the whole acces- 
sible space. The counterion density at the cell boundary 
will increase for decreasing (j)p getting close to the aver- 
age density in the limit (j^p ~* 0. We remark, however, 
that the non-monotonicity occurs at high polyion packing 
fractions of order (ftp « 0.05 — 0.2 where the approxima- 
tion of a spherical cell becomes questionable. 

The PBYM for a fixed bare charge leads to larger 
values (long-dashed line in Figure 12) which still cor- 
rectly describe the trend and the non-monotonicity. If 
the SPM data for the smallest colloid concentration are 
taken as a benchmark, a bare charge of qp — 269e is 
necessary to reproduce the same renormalized charge 
within the Poisson-Boltzmann theory. This procedure is 
in strong analogy with interpreting an experiment where 
the charge is a fit parameter to describe the structural 
data. Starting from this bare charge and changing the 
colloidal density, the PBYM predicts a similar trend for 
the renormalized charge (short-dashed line in Figure 12) 
but the actual numbers are different. This is consistent 
with experiments on strongly deionized colloidal sam- 
ples which were successfully interpreted usinga Poisson- 
Boltzmann renormalized colloidal charge p^j30| , p9[ |. 

VI. EFFECTS OF ADDED SALT 

Within the HSSM, the salt ions enter as charged hard 
spheres. For simplicity, we have considered a situation 
where the salt ions are monovalent and have the same 
diameter as the counterions and the solvent. Results for 
the effective interactions for a case with added salt are 
presented in Figure |l^. As expected, the salt ions provide 
an additional screening such that the forces are less re- 
pulsive than in the salt- free case (compare with Figure Q). 
The SPM reproduces the full HSSM data for intermedi- 
ate distances but there are deviations for molecular dis- 
tances. This is in contrast to the salt-free case where 
good agreement between the SPM and the HSSM was 
found even for small distances. The physical reason for 
this is that the pair potential decomposition which is the 
basic approximation of the SPM breaks down for nearly 
touching polyions as important configurations are paired 
microions squeezed between the polyions. This is a man- 
ifest many-body situation beyond the pair level. Still for 
two well-separated polyions or a single polyion, the SPM 
and the HSSM yield similar results for the counterion 
density field or the colloidal forces. 

As in the salt-free case, we have tested the PBYM. 
The solvent-renormalized charge q* as obtained from the 
SPM is plotted as a solid line in Figure |l^ versus salt 
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concentration. It decreases for increasing salt concentra- 
tion. The PBYM for the same bare charge of Qp = 280e 
yields the same trend as obtained in earlier investigations 
Esl, see the long-dashed line in Figure 14. If scaled by 




FIG. 13. Same as in Figure y but now for added monovalent 
salt, Cs = O.O22M0I/I. 
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FIG. 14. Renormalized charge versus decadic logarithm 
logjQ (cs/co) of the salt concentration where co = lMol/1 is 
a reference salt concentration. The parameters are the same 
as in Figure |ll] but now (j>p = 8 x 10""^ and CTc = SA. Solid 
line; SPM; long dashed line: PBYM; dashed line: PBYM for 
a fixed bare charge of qp = 210e; dot-dashed line: PBYM for 
a fixed bare charge of = 130e. 

using the SPM data for the salt-free case as a bench- 
mark (short-dashed line in Figure 14), the trend obtained 



in the PBYM is almost the same as that in the SPM. 
This explains the success of fitting experimental data 
||,||j30[|2|] by using the PBYM for real colloidal sam- 
ples which typically contain a lot of added salt. If the 
SPM data for a high concentration of added salt is used 
as benchmark point, the PBYM predicts much smaller 
renormalized charge upon deionization (dot-dashed line 
in Figure 14). 



VII. COMMENTS ON OTHER MECHANISMS 
FOR POLYION-POLYION ATTRACTION 

We finally comment on two other physical mechanism 
for mutual attraction between like-charge colloids. The 
first is the counterion depletion mechanism which was 
found within the PM in salt-free colloidal suspensions 
with strong Coulomb coupling (as realized by a small 
dielectric constant) We have redone the simulation 
using the same parameters as in Rcf. |^ but now with 
added salt. The depletion attraction is reduced but still 
present, see Figure As put forward in Ref. ||], the 
range of the attraction is comparable to 



47r / a/Sctp 



(10) 



which is a typical counterion distance corresponding to 
the spacing of a triangular lattice on the spherical col- 
loidal surface. This length is also included in Figure 13. 




FIG. 15. Reduced distance-resolved force _F(r)o-p/fc_BT ver- 
sus reduced distance r jup. The counterions are divalent and 
the size asymmetry is CTp : CTc : (Ts = 33 : 1 : 1. The further 
parameters are \qplqc\ = 16, e = 5 and r^p — 1.6 x 10"^. Solid 
line: PM without salt; long-dashed line: PM with monovalent 
salt at concentration Cs = 0.003Mol/l. The range a/(Tp=0.476 
of the depletion attraction is also indicated. 
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Unfortunately, the polyion radius used in Ref. H] is too 
large to allow for a reasonable number of solvent particles 
in the solvent bath. Therefore we have slightly reduced 
the polyion size such that the PM yields the same coun- 
terion depletion-mediated attraction. Results based on 
the HSSM and PM are collected in Figure m. As can be 
deduced form this figure, the depletion-mediated attrac- 
tion is stable with respect to an explicitly added solvent. 
It is further stable but reduced with respect to added 
salt. However, an added solvent reduces the attraction 
a bit. The physical reason for that is that the solvent 
will prefer to stay in the counterion-free space near the 
colloidal surfaces such that the solvent depletion cancels 
part of the counterion depletion force. 
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FIG. 16. Reduced distance-resolved force F{r)ap/kBT 
versus reduced distance r/ap for divalent counterions and 
Gp : Gc Us ~ W : 1 : 1. The further parameters are 
\qp/qc\ = 16, e = 20 and 0p = 5.8 x 10"^ Solid line: HSSM 
without salt; long dashed line- PM result with added mono- 
valent salt at concentration Cs — 2.74 x lO^^Mol/l; dashed 
line: PM without salt; dot-dashed line: pure electrostatic in- 
teraction between a pair of -|-2e and -2e ions. The depletion 
range a/(Tp=0.476 is also shown. The long range tails of the 
forces are compared in the inset. The arrow there indicates 
the range a. 

Second, we comment on the mechanism of metastable 
oppositely-ionized colloids 
leading to long-ranged Coulomb attraction as found in 
recent salt-free PM simulations by Messina et al 0| . We 
have confirmed and reproduced this effect in our simula- 
tions for the parameters of Figure 16. We found, however, 
that the opposite ionization of the colloids will be imme- 
diately suppressed if salt is added. In Figure |l6|, it is 
shown that for large distances the salt-free PM data are 
close to the long-ranged Coulomb force for an ionization 
degree by one counterion (compare the dot-dashed and 
the short-dashed line in the inset of Figure 16). Once 



salt is added, however, the additional microscopic ions 
will be attracted directly towards the ionized polyions 
and the long-ranged Coulomb attraction disappears (see 
the long-dashed line in the inset of Figure 16). Hence, 
the mechanism of attraction due to metastable ionized 
states is not stable with respect to added salt. We also 
remark that metastable ionized states will disappear for 
separations shorter than the characteristic depletion zone 
length a. For such close configurations, the mutual at- 
tractive Coulomb correlations in the counterion cloud 
around both polyions will lead to a symmetric shearing 
of counterions by the two neighboring polyions. For such 
small separations, counterion depletion is responsible for 
the attraction. As a function of distance, the total force 
is nonmonotonic. For small distances it is dominated 
by counterion depletion which decays off rapidly on the 
scale a, while for larger distances the electrostatic result- 
ing from the metastable oppositely-ionized colloids leads 
to a long-ranged attraction. 



VIII. CONCLUSIONS 

In conclusion, based on simulations of an model which 
contains the granularity of the solvent explicitly, we 
have shown that hydration forces profoundly influence 
the colloidal interaction. For divalent counterions, there 
is solvent-induced attraction which is not contained in 
the traditional primitive model but can be encaptured 
within a solvent-averaged primitive model (SPM). For 
monovalent counterions, the forces can be described by 
a solvent-induced charge renormalization. This picture 
is in agreement with experiments on strongly deionized 
samples where a Yukawa picture can be employed pro- 
vided the colloidal charge is renormalized towards a value 
smaller than the bare charge |^^. The trends of the 
renormalized charge upon increasing the salt concentra- 
tion is similar in the Poisson-Boltzmann cell model and 
the SPM which explains why the experimental data could 
be well-described by using a Yukawa interaction with a 
Poisson-Boltzmann- renormalized charge ||2^ , |30| , ^ . Still, 
quantitatively, there are differences between the renor- 
malized charges of the Poisson-Boltzmann cell model and 
the SPM. 

Future research should focus on the role of the perma- 
nent dipole moment in a polar solvent as modelled by 
dipolar hard spheres |55 5^ or a Stockmayer liquid [ p7| . 
Also more work has to be done to explore the role of 
charge regulation and chemisorption of counterions near 
the colloidal surfaces. Furthermore the dielectric discon- 
tinuity at the colloidal surfaces resulting in image charge 
effects has to be explored in more detail. For all these cir- 
cumstances the concept of a renormalized polyion charge 
resulting in a Yukawa picture should be possible provided 
there are free counterions left which dominate the effec- 
tive repulsion between the colloids. 

We finally point out further possible applications of 
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our simulation technique: If used without the confin- 
ing solvent-bath shell, our approach starts "ab initio" 
and even employs the correct microscopic (molecular) dy- 
namics of the solvent. Therefore, it could also be used 
to address dynamical questions in equilibrium and non- 
equilibrium. Important examples concern the motion of 
poly- and counterions under the influence of an exter- 
nal electric field including effects as the electrophoretic 
mobili ty p6[ -^9[, ion migration ]50[ |, electro-kinetic prop- 
erties and electrolyte friction [^2|j5^ ]. Our approach 
produces both diffusive motion and hydrodynamic inter- 
actions mediated by the solvent as an output. Of course, 
one will not be able to simulate large time scale separa- 
tions between the Brownian and the structural relaxation 
time but one should try to start with moderate time 
scale separations in order to test the approximative the- 
ories of electrophoresis. 

We thank R. Roth, M.Schmidt and T. Palberg for help- 
ful remarks and the DFG for financial support. 
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